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Abstract 

We introduce a natural generalization of the forward-starting options, first discussed by 
M. Rubinstein ([H])- The main feature of the contract presented here is that the strike- 
determination time is not fixed ex-ante, but allowed to be random, usually related to the 
occurrence of some event, either of financial nature or not. We will call these options Ran¬ 
dom Time Forward Starting (RTFS). 

We show that, under an appropriate “martingale preserving” hypothesis, we can exhibit 
arbitrage free prices, which can be explicitly computed in many classical market models, at 
least under independence between the random time and the assets’ prices. Practical imple¬ 
mentations of the pricing methodologies are also provided. Finally a credit value adjustment 
formula for these OTC options is computed for the unilateral counterparty credit risk. 
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1 Introduction 

Forward-starting options are path dependent put/call financial contracts characterized by having 
a strike price expressed in terms of a pre-specified percentage a of the asset price taken at a 
contractually fixed intermediate date u G {t,T], T being the option maturity. The time u is 
known as strike-determination time. The payoff of a forward starting call is therefore 

(St - (1) 

These products represent the fundamental component of the so-called cliquets (see m), which 
are indeed equivalent to a series of forward starting at-the-money options, activated along a 
sequence of intermediate dates, upon payment of an initial premium. Cliquets are often employed 
to buy protection against downside risk, though preserving an upside potential, for instance in 
pension plans in order to hedge the guarantees attached to embedded equity linked products. 

Wilmott in [29] showed that these products are particularly sensitive to the model that one 
chooses for the dynamics of the underlying’s price. 
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In this paper we study a generalization of forward starting options allowing for random 
strike-determination times. That is to say, we hypothesize that we are making the forward 
starting option worth, if some future random event happens during the life of the contract with 
positive probability. The value of the asset at this random time is then compared at maturity 
with the final value of the underlying, so the call’s payoff ([I]) is replaced by 

{St - aSr^T)^, (2) 

for some appropriate random time r, where we denoted by x A y = min(x,y). The pricing 
methodology we will present is an extension of Rubinstein’s approach (see [28jl. and it can be 
applied to several existing option models. 

It is easy to imagine a promising application of this type of contracts in several situations, 
as we show with the following examples. 

• A first example with a speculative goal, could be represented by a call written (with 
a = 1) on the Manchester United stock (NYSE:MANU), starting the beginning day of 
the Premiere League (PL) soccer tournament and ending the closing day of the tournament. 
The call is triggered by the first time, r, the United will eventually lead the tournament 
with at least three points more than all other teams. If this event happens, then the United 
will increase its chances to take part to the next Champions League (CL) tournament, an 
event which generates substantial additional cash-flows. If, at the end of PL tournament, 
the United will get the right to play CL, the share value St of MANU at that date will 
probably be larger than Sr and an investor in the call contract could pocket the difference. 

• Another possible context we may imagine is when a firm expects an adverse market scenario 
that might cause a downgrading of the firm’s rating in the near future. So the firm’s owner 
offers a contract of this type to a smart manager (hired in the past), as an incentive to 
convince him/her to remain within the firm and to work to improve the firm’s position 
before maturity. As an example, the firm’s owner offers this contract to the manager on 
January, 1st 2015 with maturity T = 1 year (with a = 1). If the event {r < T} happens, 
for instance on June, 6th 2015, then the manager has time until maturity to make the 
firm’s value rise again above Sr and pocket the payoff, otherwise no reward is gained. The 
manager is therefore interested in assessing the correct current value of the incentive. 

In this context, it is also interesting to consider the counterparty credit risk aspects of 
this type of contract from the point of view of the manager: if the firm (the counterparty) 
does not recover from the adverse market scenario, it might default a short time after the 
downgrading event or default could even happen before downgrading. This leads to the 
problem of computing unilateral Credit Value Adjustment (CVA) for this kind of option, 
since the manager needs to evaluate this further risk of the contract in order to assess 
correctly the value of the firm’s offer. 

• Lastly, it is also possible to consider contracts depending on more than one single random 
event: for example a cliquet call option could be triggered by the first time the underlying 
asset reaches a pre-specified upper barrier Ki until the asset reaches, if before maturity T, 
a new barrier level K 2 > Ki , at that time a second call option with strike K 2 and expiring 
at maturity gets activated. In general, depending on the characteristics of the financial 
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product, one can imagine to construct cliquet type options along some increasing sequence 
of random times. 

The literature on classical forward starting options is wide in terms of underlying’s dynamics 
model and evaluation techniques. We shortly review the main contributions, since we will use 
these results as starting point in our analysis 

Rubinstein (1991) (see [28]) was the first to consider the pricing problem of a forward-start 
option in the Black and Scholes model. 

Later, valuation of forward starting options in stochastic volatility models was addressed by 
several authors. Lucic (2003) (see [23]) wrote an evaluation formula by employing the asset’s 
price at determination time as numeraire, transforming those options in plain vanilla ones, giving 
a closed formula for the Heston’s model. Kruse and Nogel (2005) (see [25]), followed the same 
approach, but they developed a pricing formula exploiting the distributional properties of the 
Heston’s model; they remarked the importance of such an extension (w.r.t. the Black-Scholes 
environment) since forward starting options are particularly sensitive to volatility changes. 

Further, Guo and Hung (2008) (see [19]) derived quasi-analytic pricing formulas in general 
stochastic volatility/stochastic interest rates models, presenting jumps in both the asset’s price 
and the volatility. This was important, since, in contrast to the plain vanilla case, the value of a 
forward-starting option may not always increase with maturity. It depends on the current term 
structure of interest rates. In the same context, Ahlip and Rutkowski (2009) (see [I]) chose a 
Heston’s stochastic volatility model and a CIR dynamics for the interest rates, with volatility 
and short rate both correlated with the dynamics of the stock returns. The main result is an 
analytic formula for the call case, derived by using the probabilistic approach combined with the 
Fourier inversion technique, as developed in Carr and Madan (1999) (see [IQ]). Van Haastrecht 
and Pelsser (2011), ([20]), exploiting Fourier inversion of characteristic functions, developed a 
quantitative analysis of the pricing of forward starting options under stochastic volatility and 
stochastic interest rates (with Ornstein-Uhlenbeck dynamics), confirming that these not only 
depend on future smiles, but are also very sensitive to model specihcations such as volatility, 
interest rate and correlation movements, concluding that it is of crucial importance to take all 
these factors into account for a correct valuation and risk management of these securities. More 
recently, Ramponi (2012) (see [27]) extended this Fourier analysis to regime switching models 
for the asset price. 

Finally, in 2006 Chicago Board Options Exchange introduced options on its implied volatility 
index (VIX), enhancing the interest in forward starting options to forecast future volatility. 
Albanese, Lo and Mijatovic ([2]) proposed a spectral method to evaluate numerically the joint 
distribution between the stock price and its realized variance, which gave a way of pricing 
consistently European, general accrued variance payoffs, forward-starting and VIX options. 

To evaluate our generalization of FS options to random determination time, we exploit the 
techniques developed in a number of papers (refer to [5] and to [6] for a general exposition) to 
deal with default time in credit derivatives. The main tool we employ is the so called assumption 
(H), that guarantees that all martingales remain as such when adding the information generated 
by the default time. In other words, this means that the no-arbitrage structure is not altered 
by the introduction of the random time. This leads to the so called key lemma, that allows to 
switch from the global information (generated by the asset’s price and the random time) to the 
information generated only by the observable prices. 
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Following the intensity based approach, we give a general expression of the pricing formula 
either in conditions of independence or not. Under independence between the asset’s price and 
the random time, we specialize and explicit our prices in several computable models (Black 
& Scholes, Heston, Variance Gamma) and we suggest an affine dependency structure for the 
intensity when the independence assumption is dropped, showing it is still possible to recover a 
partially easily computable formula. 

Finally, keeping in mind the second example presented above, we focus on computing the 
Credit Value Adjustment, due to counterparty risk, required by these products, when in presence 
of a further default time (for a general reference see Brigo and Chourdakis (2009), [7]). 

The paper is organized as follows. Section [2] gives the set up, in Section[3]we present a general 
approach to the evaluation of this kind of product, in Section |4] we develop some explicitly 
computable cases under independence. Section [5] provides some numerical implementation of 
the derived pricing formulas, finally Section 0 treats briefly the counterparty risk implied by this 
product, when a default time for the issuer might occur. 


2 The set up 

We consider a time interval [0,T], with T denoting the maturity of the contracts. We assume 
that the market is made up of a single asset and a bond living in a probability space {Q,T,P), 
with a filtration {T( 5 }s>o, made complete and right continuous to satisfy the usual hypotheses 
(see [26]). Moreover, we assume we are in absence of arbitrage opportunity and that P is a risk 
neutral measure selected by some criterion. 

Before formalizing mathematically the new contract previously described, we need to recall 
briefly the evaluation of classical forward-start options. 

Denoting by B[t, T) the discounting process, the price of a forward-start call with percentage 
a and strike-determination time u > t is given by 

c(t,u,r) =E^[P(t,r)(5r-aS,)+|Pt], 0<t<T. (3) 

In market models where the price of a plain vanilla call option is represented by a deterministic 
homogeneous function of degree 1 in two spatial variables {x,y), that is 

call{t, 7 x, 7 y, T) = 'ycall{t, x,y,T), V 7 , x, y > 0 


applied to the given pair the price of a classical (call) forward-start option is 

cit,u,T) = E^[B{t,u)EP[B{u,T){ST-aSu)+\Bu]\J^t] = ^^[B{t,u)call{u,Su,Sua,T)\Pt] 


= E^[B{t,u)caU{u,'yl,'ya,T) {Pt] = P‘^[B{t,u){'ycaU{u,l,a,T)) 

'y=Su 

= E^[B{t,u)SuCall{u,l,a,T)\iFt]- 


'y=Su 


\Tt] 


(4) 


Scale-invariant models (see 0) certainly hold this property. Furthermore, if the call price is 
a deterministic function of the current price of the underlying and it does not include any 
additional stochastic factor, we may conclude that the price of the forward-start option verifies 

E^[B{t,u)SuCall{u,l,a,T)\Pt] = E^[B{t,u)Su\B't]caU{u,l,a,T) = St caU{u,l,a,T), (5) 
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where in the last equality we used the martingale property of the discounted asset price. Notice 
that in such a case the portfolio having call{u, 1, a, T) units of the stock ‘replicates” the t-value 
of the forward-start option. 

Sometimes an analogous result may be achieved also in presence of stochastic factors via 
change of numeraire arriving at a formula of the type 

¥.^[B{t,u)SuCall{u,l,a,T)\Ft] = StE^[call{u,l,a,T)\Ft\ ( 6 ) 

for some proper equivalent measure P. 

We want to generalize these contracts allowing the strike-determination time to be random, 
thus we need to consider a market model (and hence a probability space) where the asset price 
S lives together with a random time r. We therefore assume that there exists a fi—algebra Q 
containing F and a filtration { 0 s}s>o ^ G-, satisfying the usual hypotheses, that makes St and 
Hi = l{T<t} adapted processes (that is to say that r is stopping time w.r.t. { 0 s}s>o ) and that 
we may extend the probability P to a new probability measure Q on the measure space (O, Q). 
Hence, referring to the previous notation, we have Ft P Gt, for all t > 0 and from now on we 
will assume Gt = B'tV Pt, where Pt = cr{Hs ■ s < t). We recall that the filtration stopped at 
the stopping time r will be denoted Gt and it means 

Gt — G 0 H n {t < t} G Gti Vt P 0}. 

Besides we introduce the classical and fundamental hypothesis 
(H) Every Pf— martingale remains a martingale. 

This assumption, known as the H-hypothesis, is widely used in the credit risk literature (see 
e.g. [IT] and references therein). With this hypothesis, we may conclude that Su remains in 
general a semimartingale and B{t,u)Su a martingale under the extended probability Q. So 
we can replace the arbitrage free pricing formula (|3|) by the more general expression 

cit,T) = E«[P(t,r)(5T - aSrAT)^\Gt]. (7) 

From now on, this product will be called a Random Time Forward Starting option and 
we will refer to it as RTFS options. Notice that the payoff remains ^'r~™easurable. In what 
follows we treat only the call case as, under hypothesis (H), a natural extension of the put-call 
parity holds for these products 

F«[P(t,T)(ST - aSrAT)^\Gt] - BQ[B{t,T){aSrAT - ST)^\Gt] = St - aSrAtB{t,T At). 

Lastly we remark that for a = 1, these options are worth something only if the random 
occurrence happens before maturity, in a sense we could say they are triggered by the random 
time. 


3 Pricing formulas 

In this section we want to understand how to make formula (|7|) more explicitly computable. 
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As a general consideration, let us notice that 


cit,T) = EQ[i?(t,r)(5T-a5,Ar)+(l{o<r<r} + l{r>T})|ai] 

= E<^[Bit,T){ST - aSr) + l{ 0 <r<T}\Gt] + (1 - a)E^[B{t,T)STlir>T}\Gt] 

= E<^[B{t,T){ST - aSr)^\gt]l{ 0 <r<t} + E<^[Bit,T){ST - aSr)^l^t^r<T}\Gt] 

+ {l-a)Stl^t^,y-{l-a)EQ[B{t,T)STlit<r<T})\Gt]. ( 8 ) 

having written as 1 — l{o<r<r}- 

For the first term there is not much to say, there t < t and it behaves as a call price with 
strike price aSr, which is completely known at time t. The third term represents the guaranteed 
payoff of this contract, so we have to study the remaining terms. Hypothesis (H) is going to 
help us for the last one. We know that B{t,u)Su is a martingale under Q and that the event 
{t < T < T} G Qr (since for any t < u < T, {t < t < T} Cl {t < u} = {t < t < u} G Qu), so 
conditioning internally w.r.t. Qr, we obtain 

E^[B{t,T)STl{t<r<T})\Gt] = E^[E^{B{t,T)STl{t<r<T}\Gr)\Gt] 

= EQ[E<^{B{t,T)ST\Gr)l{t<r<T}\Gt] = E^ {B{t,T)Srl{t<r<T}\Gt\, 

where we used the optional sampling theorem, assuming enough integrability of the asset price 
process, as r < T almost surely. 

The second term may be rewritten as 

E^[B{t,T){ST - aSr)^l{t<r<T})\Gt] = EQ[EQ{B{t,T){ST - aSr) + l{t<r<T}\Gr)\Gt] 

= E^[B{t,T)E^{B{T,T){ST - aSr)+\Grn{t<r<T}\Gt] = E^[B {t, t)c{t, t,T)1 {t<T<T}\Gt\, 

with c defined w.r.t. Qt as in ([3]), where we extended the definition to include also u = t. 

To proceed with our evaluation we have two possibilities. Either Qt = he. r is an 
stopping time, or not. In the first case, some computable situations may occur. For instance, 
if T = inf{ti > 0 : Su > H} is a hitting time (choosing the barrier level H > Sq), we have 

c(t,T) = E^[B{t,T){ST - aH)+\Bt]l{o<r<t} + ^'^[B{t,T){ST - aH')+lp<,<T}|-Ti] 

+ (1 - a)Stl{t<r} - (1 - a)E<^[B{t,T)STl{t<r<T})\Bt\. 

Summarizing 

c(t, T) = call price(t, T, St, aH")l{o<T<t} + barrier price“^/*”'(t, T, St, aH, H) (9) 
- (1 - a)EQ[B{t,T)STl{t<r<T}m + (1 - cx)Stl{t<r}, 

where “barrier price” denotes the price of a barrier call option with strike price aH which is 
activated as soon as the threshold H is reached by the asset price. If this event does not happen 
before maturity T, at maturity the contract is not going to pay zero to the holder as in standard 
barrier contracts, but (1 — a)-percent of the terminal asset price St, since the first three pieces 
of the formula are zero. Therefore, if r > T, this random starting forward option differs from a 
plain vanilla barrier option by a specific non negative value, never exceeding (1 — a)H. 
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Other situations may be more or less complex to evaluate, depending upon the definition of 
the stopping time r. 

In general, since the observable is the asset’s price process, it would be interesting to be able 
to write the pricing formula in terms of J-t, rather than Qt = J-t y T-Lt- For that we have the 
following key Lemma (see m, El)- 


Lemma 3.1 For any integrable Q—measurable r.v. Y, the following equality holds 


eq 


Mr>t}Y\gt =Q{r>t\gt)- 


EQ 




Q(r > t\J^t) 

Applying this lemma to the second and fourth term of dS]), respectively with 
Y = B{t, T)c(r, r, r)l{t<^<r} and Y = B{t, T)Srl{t<r<T} 
and remembering that 1 — Ht = l{T>t} is measurable, we obtain 

B^[B{t,T)c{T,T,T)l{t^^<T}\gt] = l{rM}- ^ - 

E [B[t,T)Srl{t<T<T}\gt\ — l{r>t} Q{t > tlBt) 


( 10 ) 


( 11 ) 

( 12 ) 


The above may be rewritten following the hazard process approach. Let us denote the conditional 
distribution of the default time r given Bt as 


Ft = Q{t <t\Ft), yt>o 

and let us notice that for u>t, Q{t < u|Ti) = E‘3((5(r < u\F'u)\iFt) = E‘3(F„|T'^). 

In order to apply the so called intensity based approach, we need to assume that Ft{u}) < 1 
for all t > 0 (which automatically excludes that gt = Bt) to well define the so called risk or 
hazard process 

Ti :=-ln(l - Et) => E* = 1 - e“’"* Vt>0, Fq = 0. (13) 


With this notation, we rewrite dill) and ()12p as 

E«[E(Lr)c(r,r,r)l{i<,<T}|et] = (1 - Ht)E<^[B{t,T)c{T,T,T){HT - Ht)eY\Bt] (14) 
EQ[B{t,r)Srlit<r<T}\gt = (1 - Ht)BQ[B{t,T)Sr{HT - Ht)eY\Tt], (15) 


but by proposition 5.1.1 (ii) page 147 of [6] these conditional expectations may be written as 


BQ[B{t,T)c{T,T,T){HT - Ht)\Bt] = 

EQ[B{t,T)Sr{HT-Ht)\Bt] = 




B{t, u)c{u, u, T)dFu |Ei] 


B{t,u)SudFu\F't] 


(16) 

( 17 ) 


If we have the additional hypotheses gt = Bt V FLt = iFt ® Fit for all t and 
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T and are independent, 


(IND) 

then r is independent of Tt-, hypothesis (H) is automatically satisfied and Ft = Q{t < t) is 
deterministic, so the previous formula becomes 

B^[B{t,T)ciT,T,T){HT - Ht)\Ft] B{t,u)c{u,u,T)dFu\Ft] 

= r BQ[-EQ{B{t,T){ST - aSu)^\Fu)\Ft]dF^ = T E«[S(t, r)(ST - a5J+| 


= J c{t,u,T)dFu 

and the other is treated similarly. 

We can summarize the above results in the following 

Proposition 3.1 With the above notation, we have 

(a) Under hypothesis (H) the price of a RTFS call option is given by 

c{t,T) = BQ[B{t,T){ST-aSr)^\gt]l{o<r<t} + {l-a)Stl{r>t} 


(18) 


+ 1 


{T>t 



B{t,u)c{u,u,T)e ^^'^dTu\Ft] 

- {1 - a)l{r^t}F^[j\{t,u)Sue-^^--^^'>dT^\Ft] 


(b) Under hypothesis {IND) the price of a RTFS call option is given by 
c{t,T) = F^[B{t,T){ST - aSr)^\Gt]l{o<r<t} 


+ 1 


{T>t} 


rT 

j c{t, u, + (1 - 


Proof: Part (a) is only the assembling of the various pieces previously presented. As for part 
(b), because of independence, P is deterministic and it gets pulled out of the expectation. Since 
Q is a martingale measure, the integrands of the last two terms verify the martingale property 
w.r.t. Ft under Q. □ 

Remark 3.1 If we assume that P^ is absolutely continuous w.r.t. the Lebesgue measure, that is 
to say 


L 


r+ = / Nds 


Vt > 0 


(19) 


for some Ft—adapted and non negative process Xt called intensity (see m), then also 

Ft = l- = 1 - e" ^0 (20) 

ft 

is absolutely continuous and if we denote by ft its density, necessarily we have Xt = -- — 

1 — It 

In this case we have 
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(a) Under hypothesis {H) the price of a RTFS call option is given by 

c{t,T) = B^[B{t,T){ST-aSr)+\gt]l{o<r<t} + {l-a)Stl{r>t} 
+ B{t,u)c{u,u,T)X^e-^t^^^‘^^du\Ft] 

- {l-a)l{r^t}^^[j\{t,u)SuXue-^t^^^^^du\Ft] 


(21) 


(b) Under hypothesis {IND) the price of a RTFS call option is given by 
c{t,T) = BQ[B{t,T){ST-aSr)^\gt]l{o<r<t} 


( 22 ) 


+ 1 


{T>t} 


c{t,u,T)Xue + {I — a)Ste 


ds 


The above formulas, either in presence of independence or not, rely on the knowledge of the 
model and of the conditional distribution of the default time. Hence, deriving an explicit or 
implementable formula for the pricing will depend heavily on the modelling choices we are going 
to make. 

Of course, the independent case is much easier than the other and we focus on it in the next 
section, to show that there are some interesting treatable cases. We will always take a = 1. 

Remark 3.2 (Affine hazard rates) 

When we no longer assume independence ofr and S, we need to resort to formula /flTl] and a 
reasonable choice is to make the hazard rate depend also on the price process. As a first attempt 
one may try an affine model, so we decide that 

Xu = a{u)Su + h{u)Zu, (23) 

where a,b are deterministic, positive and bounded functions, while Z is a positive gt—adapted 
process independent of Ft, 0 <t <T. 

By using the optional projection theorem (see m theorem VI. 57), we may arrive at the 
pricing formula 

c{t,T) = B^[B{t,T){ST - aSr)^\gt]l{o<r<t} 

rT 

+ / a{u)call{u, l,a,T)Ffi[B{t,u)Sle- 

rT 

+ {1 - a)E^[B{t,T)STe~ 

where all the above expectations are of the same type. If not explicitly computable, the expecta¬ 
tions can be evaluated via Monte Carlo or transform methods, once a grid of points t = uq < 
ui <■■■< UN = T in order to approximate the time integrals has been fixed. More precisely, 
one applies the Monte Carlo approach to the inner expectations by generating M independent 
paths of the processes S and Z at the grid points. 
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4 Some computable models 


We now show some examples where the above formulas for the RTFS options can be explicitly 
computed, of course the first case to analyze is the Black and Scholes model. 


Pricing in the Black & Scholes market 

In the classical BS market with constant risk-free rate r and volatility a, the price of a 
plain vanilla European call with maturity T and strike K will be denoted by call^^{t, St, K, T), 
while the price of of a forward-start option with maturity T, strike-determination time u and 
percentage a = 1 (see [28]) is simply 

{t,u,T) = Stcall^^ {u,l,l,cr,r,T), u>t (24) 

where 

call^^it, St,K- a, r, T) = StM{di) - KAf{d 2 ) 


di = 


log(V^) + (r + uV 2 )(r-t) 


d2 = 


log{St/K) + {r-ay2)iT-t) 


ay/T — t ay/T — t ’ 

hence in formula (p^]) we have di = cty/T — u, i = 1,2, ci = {^ + f), C 2 = ci — a. Let us 
notice that cl > r is always true. 

Henceforth, under (IND), the price of a (RTFS) option is given by 
c{t, T) = l{o<r<t}co//^‘^(t, St, Sr, a, r, T) l{r>t}Stj call^^{u, 1,1, a, r, T)Xue~ (25) 

The valuation of ([25]) can be found in closed form in some simple cases. If r ~ exp(A), A > 0 
under Q, then for the second term we obtain 

rT 


^{T>t}St call^^{u, 1,1, a, r, r)A«e = l{r>t}St Ai{t, T) - A 2 {t, T) 


(26) 


where 
All 


(t,T) = Ae^*^ M{ciVT^)e-^^du, A2{t,T) = Ae^^^ 


e-^('^-“)W(c 2 VT^)e-^“du (27) 


and it is clear that Ai > A 2 for all values of the parameters. 

The above integrals specialize depending on the choice of the parameters. If r 7 ^ A we have 
three cases. 

If Cl > 2A then these integrals can be explicitly computed exploiting integration by parts 
and the Gaussian density, obtaining (keeping in mind that c| — 2(A — r) = cf — 2A) 

Ai{t,T) =M{ciV^t)-^^==M{yJ{cl-2\){T-t)) + ^—{^^=-l) (28) 


^ci-2\ 


Vci-2A 


A2{t,T) = 


+ 


A 


A—r 

e-KT-t) 


p p — \{T —t) j - 

(^-*)Ar(c2V^) - ^^-^ M{yJ{cl-2X){T-t)) 

V Cl —2X 


(29) 


C2 




-1) 
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When cf = 2A (hence cf = 2(A — r)), then the above formulas reduce to 


A(t,T) = N{c^VT^t) - 

^ V 2vr 




A 


A — r - 


2 V 2vr 


(30) 

(31) 


and for cf < 2A we lose the Gaussian integral and we arrive at 


Ai{t,T) = MiciVr^t) 


_ p-A(r-i) 


Cl 


1 

L2 + y2A^ 

1 


d 


^{2X-ci)(T-t) 4 


zdz 


A2{t,T) = - -|e-’'('^-*W(c2V^)-e-^(^-*) 

A —r t 


+ 


C2 


V(2A-c2)(T-i) 4 


^2A —cf Jo 




dz 


(32) 
}. (33) 


If instead condition A = r holds, we arrive in any case to an explicit formula with ^i(f, T) the 
same as in (f25|) and A 2 {t,T) given by 


A2{t,T) = Ae-^(^-*) [(T - t)M{c2VT^t) + 


^ _ 1 _ 1)1 , ( 34 ) 


27rc2 


2^1 


Remark 4.1 (The Merton market model) 

The above results can be extended with moderate computational effort to the Merton jump- 
diffusion market model (see e.g. m) 

1 Q 

= (r — vK)dt + adWt + (e*^ — l)dA^t, 
h>t- 

where Wt is a Brownian motion, Nt a Poisson process with jump-arrival intensity v, J ^ 
M{pL,5‘^) all independent of each other, k = E‘3[e'^] — 1 is the jump-size expectation, a is the 
volatility and r is the risk-free rate. 

Under (IND), for a = 1 and under the martingale Q, if t ^ exp(A), with A > 0, the price of 
a (RTFS) option is given by 


■ 1 "^ jpTl 

c{t,T) = l{o<r<t}ca//'^(t,5t,5,,r) + l|,>i}5t^ —(Ai,„(t,r)-Al 2 ,„(t,r)), (35) 

n =0 

where u = i^(1 + k), Ai^n{t,T) and A 2 ^n{t-,T) are given as in |I7D for appropriate sets of param¬ 
eters varying with n. The explicit formula is available from the authors, upon request. 


Pricing with Fourier Transform: the variance Gamma model 

Here we assume that the price dynamics is of the form St = , u <T, where Xu is 

a Levy process. Fourier transform is a well established computational tool for derivative pricing, 
first introduced in [21] and then further developed by several authors (e.g. see [22| or m and 
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the references therein). Different ways of appiying this technique iead to different representation 
formuias for the risk neutrai price of a piain vaniiia caii option, we choose to consider 


call{u, Sy 


1 f 


1^+00 (j)XT-Xu{-z) 

iz — ’ 


(36) 


for z € C and u = Im{z) > 1, with cpxT-Xui^) = the generaiized characteristic 

function of Xt — Xu. Because of the independence of the increments, 4>Xt-Xu does not depend 
on Su, so the modei is scaie-invariant and, appiying ([S]), we may conciude that the price of the 
standard Forward Starting option with determination time u is 


c(t, M, T) = Stcall{u, 1, a, T) 




a r+°° 

Jiu—oo 


g—iz(— log(o)+r(T-n)) 4>Xt—Xu{ (l~ 

iz — Z^ 


By inserting this Fourier representation of the Forward Starting caii price into the formuias of 
Proposition 1.1 or Remark 13.11 we get a doubie integrai representation of our RTFS caii price. 
As usuai, under hypothesis (IND), setting a = 1, we obtain 


c(t,T) = E‘^lB(t,T)(ST - Sr) + IQt]l{0<r<t} 


+ 1 


{T>t} 


-/ 

2’rA 


iZ/'+CXD 


IZ — L 






du 


dz. 


In some cases the inner integrai can be soived in ciosed form. As a matter of fact, if we choose 
a Variance Gamma modei for the price dynamics 

Xt = bYt + cWy^ , St = , 


with Wt a standard Brownian motion and Yt a Gamma process independent of Wt, with param¬ 
eters 1 and /i and io = l//iiog(l — bfi — fic^/2) (see [H]), we have 

T —u 2 

M _ -jt{T-u)ln{l-ibfiz+c^ 11^) 


4>Xt-Xu{z) = (- 


— bfj,z -)- c^fi- 


where inr/ = in|r/| + iaig^rj), with rj £ C and —vr < arg(r]) < vr, denotes the principai compiex 
iogarithm. Taking aiso (for ease of computation) a constant hazard rate A, we may conciude 


c{t,T) = EQ[B{t,T){ST - Sr)+\Gt]l{o<r<t} 

C /•it'+oo 1 


I _ ^-lrG+iz:)+jt\n{l-ibiiz+c^iiZ^)-X]{T-t) 

-oo '^z — z"^ r{l + iz) + - ibij,z + - X 


dz, (37) 


ieading to a computabie formuia. It is reasonabie that the same arguments might be extended 
to tempered stabie distributions as weii as to the Garr-Geman-Madan-Yor modei (see IS])- 


Pricing in the Heston stochastic volatility market. 
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Assuming independence, we can arrive to an explicit formula also in the Heston stochastic 
volatility model, given (under the risk neutral measure) by 


dSt = Si(r + V^dlH/) 

dat = k {0 — cjt)dt + c^/a^{pdWl + -^/l — p^dW^), 


where r is a constant interest rate. Here Tt has to be the natural filtration generated by the 
couple {S,a), which is jointly a bidimensional Markov process, hence the final pricing formula 
will depend only upon the initial values of both S and a. 

As before, we start from the expression 


CHit,T) = B^[B{t,T)iST-aSrV\gt]l{o<r<t} 

- rT ^ ^ 


ds 


(38) 


By employing the results in [25] (or |24] ). we can find an explicit formula for the forward start 
option. Indeed by virtue of the scale invariant property verified also by this model, we have 


r+oo 

CH{t,u,T,St,at) = / StcallH{u,l,a,a,T)f^^\^^{a)da, 
Jo 


(39) 


where callniu, 1, a, cr, T) is the standard price of a call option in the Heston model (for details 
see for instance |25| Lemma 2.1) and 


fcTuWt 


B(u) b{u)it+a{u) I'B(u)cr\W^-^)/^ 
A^/ 


lR/2-i{\/Hu)B{u)a) 


B{u) = 4 (^ pc) 


A{u) = R=^, 


(40) 

(41) 

(42) 


and / is a modified Bessel function of the first kind. 


Remark 4.2 The previous method can he adapted to cover also the case of affine stochastic 
interest rate market model proposed in ]20^ . where the interest rate follows a Hull and White 
model and the stochastic volatility is a Ornstein- Uhlenbeck process. 


5 Numerical results 

In this section we report the results of the numerical implementation of some of the models 
presented before: the Black & Scholes, the Variance Gamma and the Heston models, with a 
random time r ~ exp (A) with parameter A > 0. In all experiments we set a = 1, t = 0 
and we assume the independent hypothesis (IND). Analytical or closed-form prices, eventually 
requiring numerical quadrature, are compared with standard Monte Carlo approximations of the 
arbitrage free pricing formula d?]) . Sample paths of the underlying dynamics are evaluated at the 
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random time r and the corresponding payoffs are collected to estimate ([7]). In our experiments 
we simulated M = 1000000 samples/path^. 

Adaptive Gauss-Lobatto quadrature was used to approximate representation integrals (VG 
and Heston models), which required the evaluation of extended Fourier transforms involving 
multivalued functions, such as the complex logarithm. To avoid the well-known numerical 
instabilities due to a wrong formulation of the characteristic function, we used the so-called 
rotation count algorithm by Kahl and Jackel (see e.g. [23]). All routines were coded and 
implemented in MatLab©, version 8.0 (R2012b) on an Intel Core i7 2.40 GHz machine running 
under Windows 7 with 8 GB physical memory. 

For each model we report the closed-form (CF) and the MG prices for different values of 
A, fixing the other model parameters. In parenthesis we report the absolute error w.r.t. the 
CF price, followed by the corresponding length of the 95% confidence interval. Without loss of 
generality we set Sq = 100, T = 2 and r = 0. 

In the Black &: Scholes model the RTFS price is given by formulas (IZSD-dMD. Exact sim¬ 
ulations of the geometrical Brownian motion with a = 0.2 at the random time r are used to 
estimate the pricing formulas. 



CF Price 

MC C.I. Length 

A = 0.25 

3.0989 

3.0998 (8.4 X 10-^) 3.7 x 10"^ 

A = 0.75 

6.6457 

6.6512 (5.4 X 10“3) 5.5 x 10“2 

A = 1.25 

8.3710 

8.3699 (1.1 X 10“3) 6.1 x 10“^ 

A = 1.75 

9.2709 

9.2872 (1.6 X 10-2) 6.5 x lO-^ 


Table 1: Black & Scholes model for different values of A. 


The RTFS prices in the VG model are given by formula (1371) . In this example we consider the 
same parameters as in m,b = -o .1463, c = 0.1213 and /i = 0.1686. Sample paths are simulated 
as a Gamma time-changed Brownian motion with a discretization step equal to 1 /\/M. 



CF Price 

MC C.I. Length 

A = 0.25 

2.0159 

2.0163 (4.1 X 10“^) 5.5 x lO'^ 

A = 0.75 

4.3394 

4.3426 (3.2 x lO-^) 4.1 x lO'^ 

A = 1.25 

5.4801 

5.4748 (5.3 x 10-^) 3.9 x lO'^ 

A = 1.75 

6.0796 

6.0857 (6.1 X 10-3) 5.8 x 10-^ 


Table 2: Variance Gamma model for different values of A. 


In Heston’s stochastic volatility model parameters are chosen as in [25|: do = 0.09, k = 4, 
9 = 0.06, c = 0.65 and p = —0.9. We need to compute numerically formula ()38p . which 
requires evaluating the function given by formula (|39p . This function, though theoretically 
well defined for all values of tt, generates a numerical singularity in the computations, that 
we treat by approximating ch{0,u,T, Sq, ao) with SocallH{u,l,a,E{au),T). This gives some 

^Even though we are not employing them here, we remark that variance reduction techniques (control vari¬ 
ates, antithetic variates) or more efficient simulation algorithms (e.g. the Broadie and Kaya exact simulation, 
as proposed in m) could improve the MC estimates. Furthermore, alternative simulation schemes exploiting 
Proposition O and Remark O can be easily designed. 
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Option price 


B & S model Variance Gamma model Heston model 





Figure 1: Prices of RTFS (red dash-dotted lines) and FS (blue solid lines) options in the Black & Scholes, 
Variance Gamma and Heston model as a function of the expectation of the random time r. 


advantage since we explicitly have E(cru) = (Toe~'^“ + 0(1 —Comparison with Monte Carlo 
approximations is run with sample paths simulated by the basic Euler scheme with a step length 
equal to \/M. 



CF price 

MC C.I. Length 

A = 0.25 

3.4907 

3.4903 (3.7 X lO'^) 3.6 x 10"^ 

A = 0.75 

7.5290 

7.5227 (6.2 x 10“^) 5.2 x lO"^ 

A = 1.25 

9.5307 

9.5157 (1.4 X 10“2) 5.7 x lO"^ 

A = 1.75 

10.5988 

10.5920 (6.8 X 10"3) 6.0 x lO"^ 


Table 3: Heston model for different values of A. 


Lastly, we choose the random time r ~ exp(A) with varying A > — in order to compare at time 

t = 0, the RTFS price and the price of a FS option with determination time u = E ir) = HT. 
In other words we compare E((5r — SrAT)~^) with E{{St — 5'e(^))'''), see Figure ([I]). In all the 
instances prices decrease as u increases, moreover we notice that; i) for u ^ 0 (i.e. A —>■ -|-oo) 
FS and RTFS prices both tend to the price of an ATM plain vanilla call; ii) for u —>• T (i.e. 
A —)• 1/2) the FS value trivially goes to 0, while the value of a RTFS option stays positive since 
Q(r <r)^ 1-1/e «63%. 

6 Counterparty credit risk for RTFS options 

In this section we are interested in computing the credit value adjustment (CVA) for a RTFS 
option, when one of the parties is subject to a possibly additional default time. 

The interest lies in the fact that those products would be OTC and knowing that the CVA 
is computable might be convenient. 
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We first recall what the CVA is for an additive cashflow. Let us denote by n(t, T) the 
discounted value of the cashflow between time t and T; this value has to be additive, which 
means that for times t < s < u it holds 

s) + B{t, s)Il{s,u) = (43) 

All European contingent claims with integrable, —measurable payoff h, trivially verify (|43l) 
since n(t, s) = We call Net Present Value of II 

NPY{t) = E{U{t,T)\gt). (44) 

We now look into unilateral counterparty credit risk. Let us suppose that the cashflow is between 
two parties B (buyer) and the counterparty C, that is subject to default, with default time tq- 
Then, the discounted value of the cashflow has to be adjusted by subtracting from the defaultless 
case the quantity 


CVA(t) = (1 - i?)E(l{,<,^<r}i?(t,rc) NPV(rc)+|ei) (45) 

where 0 < i? < 1 represents a deterministic recovery rate (see m)- 

In the case of classical Eorward Start options with fixed strike-determination time u, it is 
possible to compute directly the CVA. We denote by {"Hp} the filtration generated by the process 
iLp := and by = T) V . Under independence between {Ft} and {^p}, taking for 

granted that B{t, u)Su remains a martingale in the enlarged filtration {^p} w.r.t. a risk neutral 
probability measure which extends the risk neutral probability P, we have 

CVA^^(t) = (1 - R)Q^{t <Tc< {B{t,T){ST - aSu)+\Ft), (46) 

where the key Lemma and the same notation as in section 1 were used, under the assumption 
that Q^{tc > t) = 1. 

To extend the previous formula to the case of the RTES options we need to consider two 
random times: the strike-determination time r and the counterparty default time tc- 

Hence we need to extend the original filtration {Ft} with both {Bt} and {"Hp}, obtaining 

Gt = Ft^Bt^ B^f 

Accordingly, we assume there exists an extension of the risk neutral probability P to a common 
probability space (H, Qt) that we still denote by Q and we set the corresponding key assumption 

(HC) Every Ft martingale remains a Qt martingale. 

Taking a = 1 for simplicity, we have that the price of a defaultable RTES call option is 

c(t,r) := E^[l{^^^T}B{t,T){ST - SrAT)^\gt]. 

So, in presence of a recovery rate R, we obtain that the CVA for this product has to be 

CVA(t,r) := {l-R)[c{t,T)-cit,T)] = il-R)BQ[l{t<rc<T}Bit,T){ST-SrAT)^\gt]. (47) 
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We may decompose the expectation as 


B^[l{t<rc<T}Bit,T){ST - Sr^T)^\gt] 

= EQ[l{i<,^<T}l{i<.<r}S(t,r)(5T - Sr)^\gt] + l{r<t}^^[l{t<rcmB{t,T)iST - 5^+1^*], 

where the last equality is justified by the fact that the first term in the first passage is equal 
to 0. In the second equality the crucial term is the first, as the second reduces to the CVA of 
a standard forward start option, for which formula (14611 applies. We may handle the first term 
under condition of independence between tc and (r, {S'*}). Denoting by = W* V Tit, we have 
gt = V and applying the key lemma to gt and we obtain 

EQm,^,^^T}l{t<r<T}B{t,T){ST - Sr)+\gt] 

,^^[Mt<rc<T}Ht<r<T}B{t,T){ST - S.)+|.F/^] 

= Q{Tc>t\gt) -—— I -• 

q{tc > 

The independence of from the remaining factors gives 

E^[l{t<rc<T}Mt<r<T}B{t,T){ST - Sr)+\gt] 

= Q{rc > - Sr)+\F^] 

and if Q{tc > t) = 1, the above reduces to 

E«[l{i<.^<T}l{i<r<r}i?(Tr)(5T-5.)+|et] = Q{t <rc< T)E^[l{t<r<T}B{t,T){ST-Sr)^\Fi^] 

which is the weighted value of a RTFS option, exactly as it happened with the standard FS 
option. Summarizing we have 

Proposition 6.1 Under hypothesis (HC) and the condition that Q{tc > t) = 1, if tc is in¬ 
dependent of (r, {Sf}), then the CVA of a defaultable RTFS call option of price c{t,T), with 
recovery rate R € [0,1], is given by CVA{t,T) = (1 — R)Q(t < tc < T)c{t,T). 

Lastly we remark that if the two random times coincide, Q(r = tq) = 1, then we may 
conclude that 

CVA{t,T) = (1 -R)c(t,r), 
which implies that the defaultable price is c{t,T) = Rc{t,T). 
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